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Blue phases are networks of disclination lines, which occur in cholesteric liquid crystals near the transition to the isotropic 
phase. They have recently been used for the new generation of fast switching liquid crystal displays. Here we study numerically 
the steady states and switching hydrodynamics of blue phase I (BPI) and blue phase II (BPII) cells subjected to an electric field. 
When the field is on, there are three regimes: for very weak fields (and strong anchoring at the boundaries) the blue phases are 
almost unaffected, for intermediate fields the disclinations twist (for BPI) and unzip (for BPII), whereas for very large voltages 
the network dissolves in the bulk of the cell. Interestingly, we find that a BPII cell can recover its original structure when the 
field is switched off, whereas a BPI cell is found to be trapped more easily into metastable configurations. The kinetic pathways 
followed during switching on and off entails dramatic reorganisation of the discli nation networks. We also discuss the effect of 
changing the director field anchoring at the boundary planes and of varying the direction of the applied field. 



1 Introduction 

In many cholesteric liquid crystals (CLC) the transition be- 
tween the regular cholesteric phase and the isotropic phase oc- 
curs through a cascade of intermediate phases known as "blue 
phases" (BP) 1 . Blue phases occur because the locally ener- 
getically favoured double twisting (i.e. local twisting in two 
directions, as opposed to the single twist of the cholesteric) is 
incompatible with the requirement of continuity 1 . This results 
in the formation of networks of defects that separate isolated 
stable regions. The network can lead to quite intricate struc- 
tures that are periodic at length scales comparable to the wave- 
length of light - this fact is at the origin of their name (indeed 
blue phases can be of many different colours). 
\ For example the BPI and BPII phases have a 3D cubic 
orientational order with lattice periods of several hundred 
nanometres, and they exhibit selective Bragg reflections in the 
range of visible light. For this reason they have interesting ap- 
plications in fast light modulators 2,3 or tunable photonic crys- 
tals^. Blue phases are in general thermodynamically stable 
over a narrow temperature range but recently their stabiliza- 
tion over a range of more than 60 K, including room temper- 
ature, have been shown 5- -. This has opened up new avenues 
in liquid crystal technologies culminating in the recent fabri- 
cation of a blue-phase display device with very fast switching 
and response times 9 . 
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Blue phases are fully 3D structures and the understanding of 
even their statics is quite a difficult task to achieve that requires 
very refined techniques and important computational efforts. 
For these reasons theoretical investigations of blue phases in 
presence of an electric field have been essentially restricted to 
the case of infinite samples. Furthermore, typically analytical 
and semi-analytical theories have focused on the case in which 
the electric field is infinitesimally smal l 1Q i n . This is in con- 
trast with the experimental characterization of the known blue 
phases in electric and magnetic fields^ - —. 

Therefore here we aim to investigate the switching dynam- 
ics of cholesteric blue phases in an electric field, by numeri- 
cally solving the hydrodynamic equations of motion for BPI 
and BPII. The choice of a realistic boundary condition is far 
from trivial in our system, and in this work we will compare 
and discuss different possibilities. We first focus on the cases 
in which the director field is fixed at the top and bottom bound- 
aries (z = and z = L respectively, where L is the sample size), 
to its stable structure in the absence of a field - this may in 
practice be achieved by pinning due to impurities or via sur- 
face "memory" effects Q We consider the switching dynamics 
both with an electric (or magnetic) field can be switched on ei- 
ther along the [0,0, 1] z direction or along the [1,1,1] diagonal 
of the unit cell (Section IIIB). Since in practice the orienta- 
tion of the director field at the boundaries can be controlled 
and varied (for example by rubbing the plates or treating them 
chemically), we then investigate the effect that different align- 

* Note that similar fixed boundaries have been used and have proved very use- 
ful to get physical insights into, for instance, permeation flows in cholester- 
ic^kml, where they lead to better comparison with experiments as opposed 
to other boundary conditions. We will discuss boundary conditions and their 
roles more in detail when commenting the results later on. 



ments conditions at the boundaries may have on the structure 
in the absence and presence of a field (Section IIIC). With re- 
spect to previous work on blue phases in electric fields—"—, 
our work considers hydrodynamic interactions, and focuses on 
switching on and off dynamics in a confined as opposed to a 
bulk system. This case should be relevant to devices which are 
typically in the jum range. 

A general picture that emerges out of this study is that the 
switching is intimately interconnected with the dynamics of 
the disclination network constituting these phases. For exam- 
ple we show that, under a small field, the disclination network 
of BPI twists and that of BPII stretches and unzips at the cen- 
ter. Above a given threshold of the voltage, in both cases the 
network breaks down. The disclinations then pin to the bound- 
aries (we restrict ourselves to the case in which the BP texture 
is anchored at the top and bottom plates). For BPII we observe 
that the breaking and pinning dynamics involve the creation of 
a disclination loop in the middle of the sample. This indeed 
allows the rotation of the disclination pattern at the boundaries 
preserving its continuity. As the field is switched off, the struc- 
ture goes back to the original state or to a cholesteric structure 
according to the steady state attained with the field on. Mea- 
surements of the time evolution of the free-energy of the sys- 
tem allow to characterize the intermediate states reached by 
the system during the switching on and off cycles. 

2 Models and methods 

Equilibrium properties of blue phases are described by a 
Landau-de Gennes free energy density written in terms of the 
tensor order parameter <2 a |3 1 i 19 i 27 . This comprises a bulk term 
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where K is the elastic constant and the pitch of the cholesteric 
liquid crystal is given by p = 2%/qo. The tensor 8 a ^§ is the 
Levi-Civita antisymmetric third-rank tensor, Ao is a constant 
and y controls the magnitude of order (it plays the role of an 
effective temperature for thermotropic liquid crystals). No- 
tice that the magnitude of order, q(r) (which lies between 
and 1), is 3/2 times the largest eigenvalue of Q. 

The anchoring of the director field on the boundary surfaces 
to a chosen director h° is ensured by adding the surface term 
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where the parameter Wo controls the anchoring strength, while 
So determines the degree of the surface order. If the surface 
order is equal to the bulk order, So should be taken equal to 
q = 1/2 (the value of q(r) in the bulk, away from disclina- 
tions). In the case of strong anchoring Wo is practically set to 
infinity, which means that the order parameter Q is fixed on 
the surfaces - this is how we imposed the anchoring in our 
simulations. 

The last contribution to the free-energy is given by 
— (e fl /127i)iS(xQ(xp^p> an d describes the interactions with an 
external electric field E, where £ a > is the dielectric 
anisotropy of the liquid crystal. The reduced potential is 
AV Z = E Z L Z , where E z is the component of the electric field 
along the x-direction and L z is the size of the cell in the same 
direction. Similar expressions can be written for the x and the 
y components. The time evolution of the system is then gov- 
erned by the equation of motion for the order parameter 18 ^ 8 , 
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where T is a collective rotational diffusion constant. The first 
term on the left-hand side of Eq. © is the material derivative 
describing the usual time dependence of a quantity advected 
by a fluid with velocity u. Since for rod-like molecules the or- 
der parameter distribution can be both rotated and stretched by 
flow gradients 28 a second term must be added to the material 
derivative: 

S(W,Q) = (§D + (D)(Q + I/3) + (Q + I/3)(§D-(D) 
2^(Q + I/3)Tr(QW) (6) 

where Tr denotes the tensorial trace, while D = (W + W r ) /2 
and CO = (W — W r )/2 are the symmetric and the anti- 
symmetric part of the velocity gradient tensor W a $ = d$u a . 
The constant £ depends on the aspect ratio of the molecules 
forming a given liquid crystal. This parameter controls 
whether the director field is flow aligning in shear flow (^ > 
2^=), creating a stable response, as opposed to flow tumbling, 

which gives an unstable response < ^=). 

The molecular field H, which provides the driving motion 
towards the minimum of the free energy, is given by 
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The fluid velocity, u, obeys the continuity equation and the 
Navier-Stokes equation, which in the incompressibile fluid 
limit reads 



p(3, + upd$)u a = 3p(IT a p) +r|3p(3 a t/p + 3pw a ). (8) 
with a stress tensor generalized to describe the flow of 



cholesteric liquid crystals, and equal to 
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where p is the fluid density and T| an isotropic viscosity. Po is 
a constant in the simulations reported here. 

Notice that, unless the flow field is zero (u = 0), the dynam- 
ics of the order parameter are not purely relaxational. Con- 
versely, the order parameter field affects the dynamics of the 
flow field through the stress tensor, which depends on the 
molecular field H and on Q. This is the back- flow coupling. 
In our simulations, the coupling to hydrodynamics via a non- 
trivial pressure tensor can be switched off, essentially by im- 
posing a constant zero velocity profile. In this way the effects 
of flow and backflow may be easily tested. 

To solve these equations we use a 3D hybrid lattice Boltz- 
mann algorithm that consists in solving Eq. © via a finite 
difference predictor corrector algorithm while the integration 
of the Navier- Stokes equation ([8]) is taken care by a standard 
Lattice Boltzmann (LB) approach. The hybrid method has 
the advantage that involves consistently smaller memory re- 
quirements with respect to the full LB counterpart^^ and its 
efficiency has been indeed already tested for binary fluids—, 
nematic devices^ and active nematic liquid crystals^ - —. 

In our simulations we fixed T = 0.3, £ = 0.7 and T| = 1.333 
as suggested by previous experience on numerical investiga- 
tion of blue phases 36 . Typical simulations reported here re- 
quired about 2000 single processor cpu hours. Most of our 
runs were parallel (MPI architecture) with 8-16 nodes. 

To locate the defects in the simulation domain, we look at 
the local behaviour of the order, q(f): namely when q(r) falls 
below a predetermined threshold, we identifiy that lattice point 
as belonging to a disclination. This prescription is very easy to 
implement and allows a determination of the disclinations that 
turned out to be as accurate as the more complicated procedure 
based on the degree of biaxiality— . For our study we have 
chosen a threshold of q = 0.19 that corresponds to the 60% 
of the largest eigenvalue of Q at the steady state. Clearly a 
variation of the threshold value generally leads to a variation 
in the thickness of the rendering of the disclinations in the 
figures. 

In our calculations the size of a cubic unit cell has been set 
equal to 32x32x32in simulation units. Note that the size 
of the unit cell which minimises the free energy is in general 
different from the cholesteric pitch - typically it is larger, and 
this size change has been referred to as "redshift" in the lit- 
eratur e 36 i 38 i 39 . In our simulations we have kept the redshift 
constant to the values suggested in Ref.—, and valid for zero 



field. We note that in principle the electric field modifies the 
optimal redshift in equilibrium—. While this is important to 
obtain an accurate phase diagram for blue phases in the bulk, 
it does not qualitatively affect either the switching dynamics 
or the behaviour of a confined blue phase in an electric field 
(we have verified this by running our simulations with no red- 
shift, both the succession of phases in Fig. 1 and the switching 
dynamics are unaffected). 

3 Results 

3.1 Steady state in an electric field 

The equilibrium structures in the absence of an electric field 
(Fig. la and Fig. Id) were obtained by relaxing the free energy 
to its minimum value by solving Eq. (|5]) numerically. Periodic 
boundary conditions were employed along the x and y axes, 
whereas on the z = 0, z = L planes we fixed the director field 
as in the bulk equilibrated blue phase configuration. Before 
the network disrupts, also the experiments in Ref. 5 suggest 
that the texture at the boundaries does not modify so that our 
assumption should be appropriate. In practice, as commonly 
assumed in the theory of permeation in cholesterics— ~— , fixed 
boundary conditions may be due to pinning of disclinations 
via surface impurities, or by surface "memory" effects (see 
below for results with different boundaries and a more thor- 
ough discussion of their roles). 

When we apply an electric field we initially consider a sys- 
tem with width L equal to the periodicity of the BP structures 
along z. Initial configurations were set according to the ap- 
proximate solutions in Ref.-. Note that the periodicity of the 
ground state is related to the blue phase pitch via a non-trivial 
constant (also known as "redshift") which we have taken to be 
equal to the one resulting from the semi-analytical minimisa- 
tion of the free energy presented in Ref. 2Q i 35 i 36 1 38 . The struc- 
ture we get then only depends on the liquid crystal chirality 
K, reduced temperature x and on £ which takes into account 
electric field effects. These are related to our parameters as 
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The static phase diagram at £ = is in a better agreement 
with experiments with respect to earlier calculations based 
on semi-analytical approximations and is given in Refs.^*2&, 
where it is shown that BPI and BPII appear in order of increas- 
ing values of k as found experimentally. 

We have chosen Aq, K and y in order to be in the appropri- 
ate region of the phase diagram. We have set Aq = 0.0034, 



K = 0.005, y = 3.775 for BPI, andA = 0.00212, K = 0.005, 
y = 2.9 for BPII. These may correspond, for instance, to a 
blue phase with lattice constant equal to 400 nm, with a rota- 
tional viscosity equal to 1 Poise and (Frank) elastic constants 
equal to — 10 pN and — 16 pN for BPI and BPII respectively, 
which are in the range of typical liquid crystalline materials. 
Accordingly, one may calculate that one space and time LB 
unit correspond to 0.0125 jam and 0.0013 jus for BPI, and to 
0.0125 jum and 0.00 17^s for BPII. Furthermore, in order to 
map electric field related quantities to real units, the best av- 
enue is to use the previously defined dimensionless number *£ 
and compare simulation and experimental values. An electric 
field of 1 Y/jum, with a dielectric anisotropy — 10, by assum- 
ing 2^ ~ 2 — 5 x 10 _6 J _1 m 3 — , would lead, for instance, to 
£ —0.01 whereas a voltage of 0.1 in Fig. 1 corresponds in our 
simulations to £ -0.1 (-0.09 for BPI and -0.13 for BPII). 

In Fig. 1 we show how the steady state, after the field is 
switched on, changes as a function of the voltage applied 
along the z direction. In the case of BPII (top row), for 
small electric field the steady state is indistinguishable from 
the starting configuration (Fig. la), while above a threshold 
(a voltage AV ~ 0.175 in simulation units with the parame- 
ters we chose, see caption to Fig. 1), the network unzips at the 
centre (Fig. lb). For a voltage larger than about 0.3, the discli- 
nations reconstruct at the surface to form an array of straight 
parallel lines (Fig. lc). In the case of BPI (bottom row), the 
zero-field structure (Fig. Id), is unstable after a first thresh- 
old (0.275) and becomes twisted in steady state with a voltage 
on (Fig. le). If the field exceeds a second threshold (about 
0.425 in the case of Fig. 1), the network again breaks and the 
disclinations pin to the boundaries (Fig. If). 

The difference in the optical signal transmitted by a liquid 
crystalline sample between cross polarizers with and without 
an applied field depends on the difference in the director field 
profile. Fig. 2 shows the director field on a cut of the sam- 
ple parallel to the xy plane at z = L/2, with different applied 
voltages. The double twist regions separating the disclination 
lines - typical of blue phases in equilibrium - are evident from 
the figure. 

3.2 Switching dynamics of blue phase devices 

We now discuss the switching on and off dynamics of a blue 
phase device - focussing on the case of strong anchoring (i.e. 
Wo » 1). Fig. 3 (snapshots to — ts) shows the response of a 
BPII cell to an electric field along the z direction. The field is 
strong enough to pin the disclinations to the surface (see the 
preceding Section and the discussion of Fig. 1). The disclina- 
tions unbind at the centre of the cell soon after the application 
of the electric field. These then join to the defect at the bound- 
aries, to reconstruct later on via the formation of a loop in the 
bulk of the sample. Such a remodelling of the defect networks 



leads to a rotation of the surface disclinations, which end up 
straight and parallel to each other. 

It is interesting to ask what happens to this highly distorted 
BP cell once the electric field is removed, especially in the 
high field regime on which we focus later on. Do we get back 
to the initial BPII configuration, or does the device get stuck 
into a different metastable configuration? To address these 
questions, we have simulated the switching off dynamics start- 
ing from the configuration shown in Fig. 3, snapshot ts. When 
the electric field is switched off to zero, the defects depin and 
meet in the bulk (snapshots %—ti). There, they join with the 
surface defects due to the anchoring to form a fully connected 
structure with the disclinations closeby in the center of the cell 
(snapshot t%). These then tilt and twist (snapshots tg — t\o), to 
finally reform the original structure (snapshot t\\). Therefore 
the BPII cell has the potential to be used as an electrically 
switchable device, provided at least that the original network 
of defects is pinned to the surface. 

Fig. 4 shows the results of an analogous electric field cycle 
for a BPI cell - note that this is the BP used in the experi- 
mentally constructed device in Ref. 7 . Our results suggest that 
switching in this disclination network both entails a qualita- 
tively different pathway and leads to a distinct steady state. 
During switching on (to —ts), the BPI disclinations transiently 
twist up, as in the steady states shown in Fig. 1. When the 
twisted defects get close enough, they recombine and once 
again pin to the boundaries. After the field is switched off 
(t6—hi), different disclinations depin and follow different dy- 
namics: some migrate back to the centre, while others re- 
join to neighbouring ones to form either arcs connected to the 
same surface, or lines spanning the whole device. In this case, 
then, the network gets stuck "en route" into a metastable state, 
which is distinct from both the zero and the high field BPI 
configurations. Despite this, it is still possible to use this BPI 
cell as a device, which switches between the configurations 
in ts and t\\ under the action of an applied field. We have 
done systematic simulations of several cycles and we have ob- 
served that there is a reversible cycle between the field induced 
state ts and the metastable state tn, but the system can not 
get back to the original BPI defect structure as at to. There- 
fore the structure is permanently deformed. It may be inter- 
esting to systematically scan the chirality-temperature phase 
space to investigate whether our finding that BPI devices do 
not switch back to their original, zero field, stable state is ro- 
bust and generic, or depends on the exact position of the start- 
ing state in the phase diagram. To answer this question, we 
have performed a simulation similar to the one described in 
Fig. 4 (i.e. field cycle for a BPI phase) but starting from a 
different point of the phase diagram. 

Fig. 5 shows the time behaviour of the free energy for the 
BPII (left) and the BPI (right) devices corresponding to the 
switching on-off dynamics of Fig. 3 and 4. In both cases the 
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Fig. 1 Steady state disclination networks for BPI and BPII under a varying field, and strong anchoring of the zero-field structure on the 
surface. BPII structures are shown in the first row, BPI are in the bottom row. The simulation domain has size 32 x 32 x 32. The range of 
stability in voltage (simulation units) of each of the networks shown in the snapshots is also indicated. Snapshots correspond to steady states 
with AV equal to (a,d), 0.3 (b), 0.4 (c); 0.375 (e), 0.5 (f). 
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Fig. 2 Director field profiles on the xy plane taken at z = L/2. Panels (a) and (b) correspond to the steady state disclination networks of 
Fig la and lb (BPII). Panels (c) and (d) correspond to Fig. Id and Fig. le (BPI). 




Fig. 3 Evolution of the disclination network of the BPII phase during the switching on (to —t$) -off (t§ —t\\) dynamics in a cell of size 
32 x 32 x 32. Strong anchoring at the z = and z = L planes are considered. The electric field, applied along the z-direction ([0, 0, 1]), is 
switched on (AV = 0.4) at to = 35000 and switched off at t$ = 150000. After a cycle, the system reforms the original defect structure, as can 
be seen at time t\\. For this reason the BPII cell can be used as a switchable device. All the values reported are in simulation units. 




Fig. 4 Evolution of the disclination network of the BPI phase during the switching on (to —ts) -off (t§ —t\\) dynamics in a cell of size 
32 x 32 x 32. Strong anchoring at the z = and z = L planes are considered. The electric field, applied along the z-direction ([0, 0, 1]), is 
switched on (AV = 0.5) at to = 35000 and switched off at t$ = 250000. The final state at time t\\ has a different defect structure respect to the 
typical network of the BPI (to). The free-energies of these two states are different (see Fig. 5) but, after several cycles, the switching is 
reversible between the states t$ and t\\. All the values reported are in simulation units. 



electric field was applied after to = 35000 iterations, and cor- 
respondingly the free energy abruptly drops. When the elec- 
tric field is switched off (t 5 = 150000 (left) and t 5 = 250000 
(right)), the free energy increases abrubtly and then relax 
down to a new (zero field) steady state. From the inset it is 
shown that whereas the BPII device is able to get back to the 
zero field state after removal of the field, the BPI cell gets 
trapped in a metastable configuration, with free energy slightly 
larger than that of the stable BPI phase. The failure of the BPI 
device to come back to its original equilibrium state seems to 
be quite generic since we observe the same phenomena for a 
starting state that belongs to a different region of the phase 
diagram (see squares in Fig. 5 right panel). 

As the blue phases are anisotropic phases, the orientation of 
the electric field matters, and we have compared the case pre- 
sented above in which the field was along the z axes to the one 
in which it is applied along the [1 1 1] diagonal of the unit cell. 
The results are shown in Figs. 6 and 7, for BPII and BPI re- 
spectively. Note that here we consider the dynamics for values 
of the field which are sufficient to break the initial networks, 
so that these can be compared with the switching dynamics 
with the field in the [001] direction. The case of a diagonal 
electric field was previously considered in Ref . 23 , although in 
that work backflow effects were not implemented and periodic 
boundary conditions were considered in all directions. 

In the case of BPII (Fig. 6), the switching on dynamics is 
now different. In particular, the network now does not unbind, 
and the central defect which joins up four disclinations is now 
stable: it moves along the direction of the field until it touches 
the surface to recombine into a pair of arcs (to — t$). Once 
more, when the field is removed the network grows from the 
surface into the bulk and reconstructs the initial configuration 
(te-hi)- 

Fig. 7 shows the switching of a BPI cell with a field along 
the [111] directions. As for the case of the field along [001], 
the disclination twists in an applied field - remarkably, the one 
along the field direction however does not (t\). When the de- 
formation becomes large enough, the network forms branch- 
ing points, which are unstable and annihilate in the bulk leav- 
ing once more a state with arcs of defects close to the bound- 
aries (t2 — t$)- this time the state is not symmetric, unlike the 
analogous case with the field along [001]. The switching off 
proceeds along a similar pathway as the one shown in Fig. 4, 
even though the final state is yet different (t§ — 1\\). 

In our simulations, we can switch off hydrodynamic inter- 
actions, which is often referred to as backflow in the field 
of liquid crystal dynamics, and therefore unambiguously pin- 
point the effect of hydrodynamics. Qualitatively, it is some- 
times the case that hydrodynamics allows the kinetics of the 
system to cross free energy barriers which would otherwise 
trap it into metastable states. Quantitatively, the relaxation 
times are usually decreased by backflow, which "speeds up" 



many processes. Fig. 8 shows the comparison between two 
switching on and off simulations, one with hydrodynamics 
and the other without - we considered the case of a BPI cell 
with field along [001]. It can be seen that hydrodynamics 
slightly accelerates the dynamics. The steady states are not af- 
fected for BPII; they are different for BPI, but are metastable 
in both cases and very close in free energy. 

3.3 Effect of boundary conditions 

We now come back to the important, and non-trivial, issue of 
the choice of the boundary conditions for our device switch- 
ing simulations. As mentioned above, fixed boundaries have 
often been used in the physics of cholesteric liquid crystals 
and blue phases-i^i^S. Indeed, in the case of permeation, 
fixed boundaries have proved to be closer to rheology experi- 
ments than, say, the case in which the director field can rotate 
freely at the boundary 18 . Still, these boundary conditions are 
very difficult to exactly enforce in experiments. In practice, 
as for cholestericsJ^ii, one may envisage that the liquid crys- 
talline texture and disclination pattern may be pinned at least 
partially by impurities, or it may be kept for some time due 
to surface "memory" effects which slow down the dynamics 
close to the boundaries. In the case of our blue phase device, 
for large electric field the regrowth of a disclination pattern af- 
ter switching off (BPII or the metastable BPI state in Fig. 4) 
is likely to be efficient due to the imprinting at the boundaries. 
To see how realistic this is, one would need experiments like 
those in Ref. 5 , but for large electric fields which destroy the 
initial blue phase at the bulk. For blue phase devices such as 
those in- the switching back occurs on the sub-ms scale: if 
surfaces respond as free boundaries, however, there should be 
a field above which regrowth is much slower, as it should then 
occur via nucleation of the stable zero-field blue phase start- 
ing from a field-induced metastable state. New experiments 
would be needed to completely settle this issue. 

On the other hand, experimentalists do have at their disposal 
a range of well developed and established techniques to accu- 
rately control the alignment of the director field at the bound- 
aries of a liquid crystalline sample or device. For instance, 
by rubbing it is possible to obtain a homogeneous anchoring, 
parallel to the boundary plane, whereas by means of chemical 
treatment the surface may be made to favour homeotropic, or 
perpendicular, anchoring. These two anchoring conditions are 
very much used in practice and, as they are very well control- 
lable, therefore it would be desirable to know what their effect 
is on the director conformation and switching of blue phase 
devices. Indeed, some simulations with relatively small sam- 
ple sizes and no hydrodynamics suggested that at least when 
the boundaries are close by homeotropic anchoring 24 , bound- 
ary conditions change the defect structure of a confined BPI 
phase. 




100000 



200000 

t 



300000 



400000 



(b) 




K=0.9,T=-1.03 
k=1.3, x=-1.98 





\ k=0.9,x=- 


1.03 




U k=i.3,t=- 


1.98 










200000 400000 600000 
t 



200000 



400000 
t 



600000 



800000 



Fig. 5 Time evolution of the free energy of the BPII (left) and BPI (right) phases during the switching on-off dynamics of Fig. 3 and 4. The 
two insets are zooms of the main figures and show the values of the free energies at the beginning and at the end of the switching on-off cycle. 
One can notice that, while for BPII the two values are essentially the same, for the BPI phase the final state has a higher values of the free 
energy (metastable state). This turns out to be true even if we start from initial states belonging to different regions of the phase diagram 
(open circles and filled squares, see legend for values of K and x). Open circles correspond to the parameters in Fig. 4; filled squares to 
A = 0.0012, K = 0.005, y= 5.085, and other parameters as above. 



To address the role of homeotropic and homogeneous an- 
choring on the structure of a BP device at zero applied field, 
we have numerically simulated the relaxation to equilibrium 
of three BPI and BPII unit cells, in a 32 x 32 x 96 simula- 
tion domain (the size was increased to minimize the effect 
of boundary layers). Periodic boundary conditions were ap- 
plied along the x and y directions, whereas (infinitely) strong 
homeotropic anchoring was enforced on the z = and z = L 
planes. Fig. 9 shows the steady states for BPII (a) and BPI 
(b) respectively. In both cases, while the bulk of the sample is 
virtually indistinguishable from an equilibrated BPI/BPII cell, 
there are strong anchoring driven deformations close to the 
boundary. In the case of BPII (panel a) these result into the 
formation of a stack of parallel disclinations confined at the 
surface and disjointed from the rest of the network. In the BPI 
cell, on the other hand, the diagonal disclinations approaching 
the boundaries bend away from these, and again more discli- 
nation lines appear at the surfaces, although this time they are 
not straight, but twisted, and the lines on the two boundaries 
are perpendicular to each other. 

Fig. 10 shows the steady state configurations found with 
homogeneous anchoring, along the z = and z = L planes. In 
this case, there are no surface defects disjointed from the bulk 
of the network, rather the network deforms and reconstructs to 
join just above the surface and avoids it. As in the homeotropic 
case, in our 3 cell domain the bulk is basically unaffected by 
the anchoring, which suggests that surface effects do not ap- 
preciably extend over one unit cell. 



We have therefore seen that changing anchoring conditions 
has a strong effect on the zero-field structure of the network 
close to the surface. What about the behaviour in an electric 
field? Rather than doing an exhaustive analysis, which is too 
costly computationally, we have chosen here to focus on one 
case, that of BPI with homeotropic anchoring. This should be 
the most interesting case, as hysteresis effects and metastabil- 
ity are more pronounced for BPI. Fig. 1 1 shows the switching 
dynamics with AV = 1 . After the field is switched on along the 
z axis, we observe that the bulk disclinations twist, as in the 
single unit cell simulation (see Fig. le), while the twisted de- 
fects appearing at the surfaces reorganize as follows: some of 
them bend and move slightly towards the bulk while the oth- 
ers remain fixed at the boundaries forming column-like struc- 
ture (to — ti). When the field is switched off (t$ —t$) the bulk 
network restructures itself, restoring the usual BPI defect dis- 
position, while at the surfaces only a partial restoration of the 
defect structure as at to occurs. Notice that for the voltage 
chosen here, AV = 1, in a device with a single unit cell, the 
disclination network breaks and does not reform (Fig. 1). In 
the current simulations, however, permanent deformations in 
the defect structure are observed at the steady state only at 
the surfaces, while the structure in the bulk is unaffected. We 
have also performed further simulations with 6 unit cells in the 
field direction (with a unit cell resolved with 16 instead of 32 
lattice points) and homogeneous and homeotropic anchoring. 
We have studied the case of a very large field which destroys 
the disclination networks. In this case the network does not re- 






Fig. 6 Evolution of the disclination network of the BPII phase during the switching on (to —ts) -off (fg —t\\) dynamics in a cell of size 
32 x 32 x 32. Strong anchoring at the z = and z = L planes are considered. Unlike the case shown in Fig. 3 the electric field is switched on 
(at to = 35000) along the diagonal ([111], pointing vertically and to the right in the picture) of the unit cell (AV/L X = AV/L y = AV/L Z = 0.2). 
The switching off occurs at t$ = 305000. As in the case in which the electric field is applied along the z-direction, here, after a cycle, the 
network disclination (t\\) is identical to zero-field configuration observed at time to. All the values reported are in simulation units. 






Fig. 7 Evolution of the disclination network of the BPI phase during the switching on (to —ts) -off (t§ —t\\) dynamics in a cell of size 
32 x 32 x 32. Strong anchoring at the z = and z = L planes are considered. Unlike the case shown in Fig. 4 the electric field is switched on 
(at t = 35000) along the diagonal ([111]) of the unit cell (AV/L X = AV/L y = AV/L Z = 0.275). The switching off occurs at t 5 = 250000. 
Similarly to the case observed in Fig. 4, after a cycle, the system can not recover the zero-field stable confioguration of the defects of the BPI. 
All the values reported are in simulation units. 
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Fig. 8 Time evolution of the free energy of the BPI phase during a switching on-off dynamics. The electric field is taken along the [001] 
direction. The stars (circles) correspond to the situation withouth (with) hydrodynamics. 



form. This may be seen as a prediction as it is possible to treat 
boundaries in liquid crystals so as to force such alignments. 
Our results also suggest that the boundary conditions are not 
expected to affect the qualitative behaviour of the devices un- 
der an electric field, at least up to intermediate values of the 
electric fields. 



4 Conclusions 

In conclusion, we have reported a numerical investigation of 
electric field induced switching of devices built starting from 
cholesteric blue phases. At steady state, with a field along 
the [001] direction (which is normal to the boundary planes), 
BPII undergoes an "unzipping transition"upon increasing the 
electric field strength; while BPI twists above a first parame- 
ter dependent electric field threshold, then breaks as the field 
exceeds a second larger critical value. It is likely that the ini- 
tial deformations of BPII (unzipping) and BPI (twisting) are 
crossovers or very smooth transitions, whereas the final defect 
disruption is a sharp transition. The dynamics of the field- 
induced disruption of the network is very rich. 

If the disclinations are fixed at the boundaries by impurities, 
our simulations predict an intricate pathway which allows the 
BPII structure to rotate its disclinations once they touch the top 
and bottom plates via the creation of a disclination loop. Once 
the field is switched off, we found that, remarkably, the liquid 
crystal can reform the original blue phase stable at zero field. 
In the case of BPI, disclinations again twist up and then pin to 
the boundary. Upon switching off the field, the disclinations 
unpin and part of the network deforms, although the zero field 
configuration is not recovered unless the applied field is below 
a critical threshold. Given the vastly different defect and di- 
rector configurations, the disclination dynamics predicted nu- 



merically should leave observable signatures in experiments. 

When the field is along the diagonal of a unit cell, i.e. along 
[111], the BPII network does not unzip, rather it smoothly 
deforms until it touches the boundary. The disclinations in 
BPI, instead, still twist and reconstruct. The different response 
to the electric fields along [001] and [111] is a result of the 
anisotropy of blue phases. 

In practice, in liquid crystal samples the orientation of the 
director field at the boundary can be controlled via rubbing or 
chemical treatment. We have therefore investigated how our 
results are affected if the director field at the surfaces of the de- 
vice is anchored and lies either perpendicularly to the surface 
or along one rubbing axis on the plane. In both cases, there 
is a boundary layer in which the disclination network recon- 
structs to either avoid the surface or form pinned disclinations 
disjoint to the rest of the network. If the device is one unit cell 
thick only, these deformations lead to new disclination net- 
works, some of which are far from the original BP networks. 
For thicker devices, on the other hand, the bulk is essentially 
indistinguistable from the thermodynamically stable phase. In 
this regime, the switching in an electric field is not very dif- 
ferent from the fixed surface boundary conditions. We also 
find that hydrodynamics, or backflow, has a relativel y minor 
effects on the switching on and off dynamics. 

We hope that our results have shown that it is now possible 
to qualitatively follow and predict the switching on and off dy- 
namics of BP devices with a variety of boundary conditions. 
The parameters we have chosen are similar to the typical liq- 
uid crystalline ones and the effects observed are expected to 
be generic and robust for different values of e.g. viscosities 
and elastic constants. It would also be interesting to study 
what happens if we relax the one elastic constant approxima- 
tion alex". Future experiments would be important to test the 
current simulation predictions and results, and to suggest fur- 
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Fig. 9 Steady state disclination networks at zero field for BPI (panel b) and BPII (panel a). The simulation domain has size 32 x 32 x 96. 
Periodic boundary conditions are considered along x and y whereas strong homeotrophic anchoring has been enforced on the planes z = and 
z = L. The effects of boundaries are limited to the size of a single domain. The bulk structure is unaffected while column-like surfaces appear 
near the boundaries. 




Fig. 10 Steady state disclination networks at zero field for BPI (panel b) and BPII (panel a). The simulation domain has size 32 x 32 x 96. 
Periodic boundary conditions are considered along x and y, and strong homogeneous anchoring along the x direction has been enforced on the 
planes z = and z = L. 




Fig. 11 Evolution of the disclination network of the BPI phase during the switching on (to — £2) -off (£3 — £5) dynamics. Homeotropic strong 
anchoring at the z = and z = L planes are considered. The electric field is switched on (AV = 1.0) at £0 = 120000 along the [001] (vertical) 
direction. The field is then switched off at £3 = 376000. 



ther technologically relevant questions. 
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